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The very high energy (VHE; E > 100 GeV) 7-ray blazar Markarian 501 has a history of extreme spectral 
variability and is an excellent laboratory for studying the physical processes within the jets of active galactic 
nuclei. A short term multi-wavelength study of Markarian 501 was coordinated in March 2009 using the Suzaku 
X-ray satellite as well as the VERITAS and MAGIC experiments to cover the critical X-ray and VHE 7-ray 
bands. The results of the quiescent-state observations with VERITAS and Suzaku are combined with public 
data from the Fermi Gamma-ray Space Telescope and compared to historical observations of the source during 
an extreme outburst, to examine the spectral variability, particularly how the spectral energy distribution varies 
with flux. 



1. Introduction 



Blazars are a sub-class of active galactic nuclei 
(AGN), and they are the dominant extra galactic pop- 
ulation in 7-rays. They have been observed to show 
rapid variability and non-thermal spectra, presenting 
a broad continuum across nearly the entire electro- 
magnetic spectrum, implying that the observed pho- 
tons originated within highly relativistic jets oriented 
very close to the observers line of sight [TS]. This 
orientation results in Doppler beaming that boosts 
the intensity and frequency of the observed jet emis- 
sion. Therefore, blazars make excellent laboratories 
for studying the physical processes within the jets of 
AGN. They were among the first sources to be de- 
tected in the VHE band and today there are 28 known 
VHE 7-ray blazars 1 . 

The spectral energy distribution (SED) of Markar- 
ian 501, a VHE 7-ray blazar, characteristically shows 
a double peaked structure. These peaks occur at keV 
and TeV energies when plotted as vFv vs v. This 
structure is common among all VHE 7-ray blazars, 
and several models have been developed to account 
for the double peaked structure. 

These models attribute the peak at keV energies to 
synchrotron radiation from relativistic electrons and 
positrons within the blazar jets. The models differ 
in accounting for the source of the TeV peak. The 
models are generally divided into two classes: lep- 
tonic and hadronic, named for their attributed source 
for the TeV peak. The leptonic models advocate in- 
verse Compton scattering to TeV energies of cither 
the synchrotron photons from within the jet or an ex- 
ternal photon field [e.g., |5J EH 03] ■ Tne hadronic 
models, however, account for the TeV emission by ttq 
or charged pion decay with subsequent synchrotron 
and/or Compton emission from decay products, or 



1 http:/ /tevcat. uchicago.edu/ 



synchrotron radiation from ultra-relativistic protons 

[e.g., mm us]. 

Observationally, Markarian 501 has been known to 
undergo both major outbursts on long time scales and 
rapid flares on short time scales, most prominently 
at keV and TeV energies. During these outbursts, 
both of the SED peaks have been observed to shift 
towards higher energies. During the most extreme 
cases, the synchrotron peak has been observed above 
100 keV. Historically, the SED has been measured in 
the VHE band only during outbursts, due to the lower 
sensitivity of the previous generations of instruments. 
This work attempts to provide state of the art short- 
term multi-wavelength measurements of the quiescent 
state of Markarian 501. These measurements are then 
compared to an extreme outburst of the source that 
was observed in 1997 |H HTj. 



2. OBSERVATIONS 
2.1. X-ray 

The Suzaku X-ray observatory, the product of a col- 
laboration between institutions in the United States 
and Japan, is an excellent tool for studying the broad- 
band spectral energy distribution (SED) of sources 
such as Markarian 501. Its multiple instruments pro- 
vide a broad energy range from 0.2 - 600 keV. Suzaku 
has two operating instruments for studying X-ray 
emission: the X-ray imaging spectrometer [XIS; [7] 
and the hard X-ray detector [HXD; [T5]. A third in- 
strument, the X-ray Spectrometer malfunctioned after 
launch and is no longer operational. 

The XIS instrument consists of four X-ray tele- 
scopes with imaging CCD cameras in their focal plane. 
Three of the CCDs are front illuminated, and one 
CCD is back illuminated to provide extended sensi- 
tivity to soft X-rays. The combined energy range of 
these CCDs is 0.2 - 12.0 keV. The HXD is a non- 
imaging instrument that expands the energy sensi- 
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tivity of Suzaku to the 10 - 600 keV band by using 
two types of detectors. Silicon PIN diodes provide 
sensitivity in the range 10 - 70 keV, and GSO scin- 
tillators placed behind the PIN diodes provide sensi- 
tivity in the range 40 - 600 keV. One limitation of 
Suzaku with respect to attaining simultaneous multi- 
wavelength observations is that it resides in a low 
earth orbit, so observations are subject to frequent 
earth occultations, limiting continuous temporal cov- 
erage. 

Observations of Markarian 501 were carried out 
with the Suzaku X-ray satellite from 2009-03-23 UT 
18:39 to 2009-03-25 UT 07:59 (sequence number 
703046010). After good time interval selection, a to- 
tal of approximately 72 ks of live time remained. The 
telescope operated with a nominal pointing that op- 
timizes the effective area of the HXD. Even with this 
pointing, the maximum useful energy for these obser- 
vations was around 40 keV. The HXD/PIN non X- 
ray background is provided by the HXD team, and 
the cosmic X-ray background was accounted for by 
estimating the PIN response to the flat emission dis- 
tribution. 

The XIS observations were requested in 1/4 win- 
dow mode to reduce the occurrence of pile up. Qual- 
ity cuts were made on the data, and XIS events were 
extracted from a circular region centered on Markar- 
ian 501. Rectangular regions far from the source were 
chosen for the background regions. The auxiliary and 
response files were created using tools developed by 
the Suzaku team. After good time interval selection, 
a total of approximately 72 ks of live time remained. 



2.2. VHE 7-ray 

VERITAS is an array of four 12 m diameter IACTs 
each focusing light on a camera consisting of 499 
PMTs. Utilizing multiple telescopes provides a stereo- 
scopic view of showers, making it possible to signifi- 
cantly reduce muon events. This substantially im- 
proves the low energy performance of the detector 
since images produced from nearby muons are other- 
wise virtually indistinguishable from images produced 
by low-energy gamma-rays. Therefore the array al- 
lows for a lower energy threshold than possible with 
a single telescope of the same size operating alone. In 
addition, stereoscopic observations facilitate a much 
improved determination of the core position of show- 
ers and thus provide improved energy resolution and 
background rejection. The increased sensitivity of 
VERITAS over previous generations of instruments 
allows for a high quality spectrum of Markarian 501 
to be collected in a relatively short time, even in the 
quiescent state. 

On 2009 March 24 2008, VERITAS took data on 
Markarian 501 while operating in wobble mode, with 
an offset angle of 0.5°. To provide background esti- 
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Figure 1: Preliminary light curves from Suzaku HXD, 
XIS and VERITAS. 



mations using this observation mode, we follow the 
reflected-region model [3] where one collection region 
is placed at the source position and others of equal 
size are placed at equal offsets from the center of the 
field of view and used for background measurements. 
VERITAS observed Markarian 501 for 2.6 hours from 
9:11 UT to 11:58 UT using a series of 20 minute runs, 
with a live time of 2.4 hours. The data were processed 
using standard VERITAS analysis cuts, and the total 
significance for this observation was 9.5 a with an av- 
erage significance per run of 3.3 a. 



2.3. HE 7-ray 

The operation of the Fermi Space Telescope allows 
for measurements to be taken in an energy range that 
has been out of reach since the Compton Gamma-ray 
Observatory de-orbited in 2000. It also has a signifi- 
cant improvement in sensitivity over previous instru- 
ments, allowing for spectral points to be measured not 
only in the flaring state, but also over the course of 
several days in the quiescent state. 

A preliminary measurement based on public data 
from the Fermi Gamma-ray Space Telescope was used 
to include a contemporaneous data point in the broad- 
band SED. There were not enough statistics from 
the two days coincident with the Suzaku observation 
to extract spectra, so the preliminary Fermi result 
shown includes data from March and April 2009 to 
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Figure 2: Preliminary VHE spectrum from VERITAS, 
including the best-fit power law. 




Figure 3: Preliminary SED with 2009 data from Suzaku 
HXD and XIS, Fermi and VERITAS (blue) as well as 
archival data from Bepposax, EGRET and Whipple 
(red). 



increase the significance of the spectral point. The 
data was analyzed following the guidelines provided 
by the Fermi team 2 . The disadvantage to the use of a 
wider time window is the inclusion of large amounts of 
data that are not simultaneous with the other obser- 
vations. At the time of writing of these proceedings, 
an updated analysis has been provided by the Fermi 
group using data from only seven days, centered on 
the Suzaku observation. This updated SED shows a 
lower flux in the Fermi energy band, and these up- 
dated results will be presented in a forthcoming pa- 
per. 

2.4. Results 

The combined data provide excellent coverage of 
the quiescent state of this extreme VHE blazar. After 
following the procedures described above and detailed 
in the ABC guide for Suzaku data analysis 3 to repro- 
cess and check the data, preliminary light curves and 
spectra were extracted. XSPEC 12 was used to fit the 
spectrum from Suzaku [2 j. A broken power law mod- 
ified by interstellar absorption fit the data well, and 
the result was unfolded and deabsorbed to derive the 
intrinsic X-ray spectrum. 

In addition, we were able to generate a light curve 
and spectrum for the VHE 7-ray data from VERITAS. 
The VHE 7-ray data was fit with a simple power law 
of the form: 



dN 
dE 



F Q • 10" 



E 



ITeV 



photons 
TeV cm 2 s 



(1) 



Finding best-fit parameters of a = 2.61 ± 0.15 and 
Fq = 5.62 ± 0.83. Uncertainties are statistical only. 

The resulting light curves are shown in Figure [l] 
In the VHE 7-ray band, the source appears to be 
in a very quiet state. For a longer term study of 
the source's behavior in the VHE 7-ray band, see [BJ. 
There was no strong variability detected by VERITAS 
or the HXD on Suzaku, and moderate variability is 
observed in the XIS light curve. 

The spectral results from the 2009 observatins dur- 
ing the quiescent state as well as archival data from 
[IT] and [I] are shown in Figure [3] 

The SED data were matched to a simple Syn- 
chrotron Sclf-Compton code written by H. Krawczyn- 
ski 4 . Here it assumes a spherical blob of radius 1 x 10 15 
cm, filled with a homogenous non-thermal electron 
population and uniform magnetic field. The electron 
population, with 7 m i„ = 1 where 7 is the Lorentz fac- 
tor of the electrons, is described by a broken power 
law. 

This simple model matches the data well, with the 
predominant change between the two flux states being 
the spectrum of injected electrons. The magnetic field 
also shifts slighty. Parameters for the models for the 
1997 flaring state and the 2009 quiescent state are 
given in Table |T] 



2 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 

3 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 4 http:/ /jelley.wustl.edu/multiwave/spectrum/?code 
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Table I SED Model Parameters 



Model Component 


1997 Model 


2009 Model 


Doppler Factor 


15 


15 


Magnetic Field (G) 


0.6 


0.95 


^max 


3.00xl0 6 


3.50xl0 5 


^Jbreak 


3.00xl0 5 


2.00xl0 4 


Index 1 (j mm - 'ybreak) 


1.9 


1.9 


Index 2 ("/break — Tmoi) 


2.6 


3.4 


normalization 


0.8 


0.8 



2.5. Discussion 

This data set provides a high quality sampling of the 
quiescent state SED of Markarian 501, and it allows 
detailed comparisons over a broad energy range to be 
made with the flaring state observed in 1997. From 
the X-ray data points alone, the synchrotron emission 
must peak above 227 keV (5.49 x 10 19 Hz) during the 
flaring state and below 0.4 keV (9.90 x 10 16 Hz) during 
the quiescent state. The SED model fits imply that 
the difference must be even greater. This means that 
with a change in flux of around one order of magnitude 
for the quiescent and flaring states used in this study, 
the synchrotron peak of the SED shifts in frequency 
by at least 3 orders of magnitude. 

The predominant change in the SED models be- 
tween the flaring and quiescent state is the spectrum 
of the electron population, implying that this may 
be the primary explanation for the dramatic shift in 
the peak frequency. The preliminary SSC model ap- 
pears to indicate that during the flaring state electrons 
are accelerated up to energies an order of magnitude 
higher than those in the quiescent state. However, 
this match is not unique due to the number of free 
parameters. An update of this study, potentially in- 
cluding more data at VHE 7-ray energies and an im- 
proved spectrum from the Fermi Gamma-ray Space 
Telescope is forthcoming. 
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